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The results of extensive Monte Carlo simulations of classical spins on the two-dimensional kagome 
lattice with only dipolar interactions are presented. In addition to revealing the six-fold degenerate 
ground state, the nature of the finite-temperature phase transition to long-range magnetic order 
is discussed. Low temperature states consisting of mixtures of degenerate ground state configura¬ 
tions separated by domain walls can be explained as a result of competing exchange-like and shape 
anisotropy-like terms in the dipolar coupling. Fluctuations between pairs of degenerate spin config¬ 
urations are found to persist well into the ordered state as the temperature is lowered until locking 
in to a low-energy state. 


I. INTRODUCTION 

The study of spin systems with short-range antifer¬ 
romagnetic exchange interactions on geometrically frus¬ 
trated lattices has revealed a remarkably wide variety 
of magnetic structures and phase transitions Exotic 
spin states resulting from frustration through the lat¬ 
tice geometry or magnetic interactions couple to other 
degrees of freedom and are associated with a range of 
disparate phenomena including anomalous hall effect)^ 
the magnetoelectric effect^ and exchange-biasing in spin 
valvesj^i^ The study of magnetic dipole-dipole interac¬ 
tions on lattice types which are known to induce frus¬ 
tration with only near-neighbor exchange coupling have 
recieved much less attention due to the relative weakness 
of magnetostatic effects. An exception is the pyrochlore 
latice and its association with spin ice phenomenal^ Re¬ 
cent attention has also been devoted to artifically tern- 
plated magnetic islands on frustrated and unfrustrated 
lattice types which are well separated and interact only 
through magnetostatic couplingiii^ 

Although the dipole-dipole interaction is much smaller 
than exchange effects, its long-range nature can lead to 
significant effects on magnetic structures in systems with 
ferromagnetic exchangei^ In cases where the fundamental 
magnetic interaction is antiferromagnetic, dipole effects 
tend to be much weaker due to smaller net magnetiza¬ 
tion. However, for low-dimensional systems and at sur¬ 
faces, there can be larger net magnetization in systems 
with antiferromagnetic exchange and dipole effects can 
be more improtant. In spin-valve structures, the antifer- 
romagnet provides exchange biasing to the adjacent thin 
ferromagnetic film, which exchange couples to itJ^ The 
impact of dipole interactions within the surface layer of 
the antiferromagnet in this case may be an important 
effect for the explanation of the exchange-bias held but 
the phenomenon of exchange bias is not well understood. 
The most popular compound for use as the antiferro¬ 
magnet in spin valves for magnetic recording is IrMna, 
which in its ordered crystalline fee phase is composed of 


ABC stacked [111] planes of magnetic Mn ions on kagome 
sites^i^di-is 

The large number of studies over the past few decades 
on the near-neighbor antiferromagnetic exchange Heisen¬ 
berg kagome model have shown that the large degener¬ 
acy associated with the 120° spin structures on corner¬ 
sharing triangles leads to the stability of so-called q = 0, 
\/3 X \/3 and other more exotic spin structures.— This 
local 120° ordering also occurs in the edge-sharing trian¬ 
gular lattice with near-neighbor exchange In the case 
of long-range dipole interactions, details of the lattice 
geometry can play a more important role. For the tri¬ 
angular lattice, the ground-state and low-temperature 
spin structures have been established to be ferromagen- 
tic whereas for the kagome lattice, this appears to not 
be the case— Spin-ice type order has been proposed 
for the kagome lattice using both Ising-like^i as well as 
anisotropic Heisenberg^^ spin models. Of particular in¬ 
terest to the present work are the Monte Carlo (MC) 
simulation studies of Tomita— on the kagome lattice us¬ 
ing a purely Heisenberg spin model with dipole interac¬ 
tions. These results suggest that alternating rows with 
six-fold ferromagnetic order, and disorder, are stabilized 
below a critical temperature of T^v — 0.43 (in dimension¬ 
less units). The ground state spin configuration was not 
studied. 

In this work, extensive Metropolis MC simulations 
are perfomed on the two-dimensional kagome Heisenberg 
spin lattice with only dipole interactions. We confirm 
a phase transition to long-range magnetic order but the 
spin configuration we find differs from that of Tomita. 
Significant fluctuations occur between coexisting six-fold 
non-ferromagnetic degenerate states over a range of tem¬ 
perature (T) before locking-in to a pair of states at lower 
T. The nature of the ground state is examined using the 
effective field method and confirmed by low-T MC simu¬ 
lation results. The impact of the range of the interaction 
on the ground state is also examined. 

The remainder of this paper is organized as follows. In 
Sec.H, the model is described and in Sec. HI the effective 
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field method (EFM) is used to reveal the six-fold degen¬ 
erate spin order in the ground state. MC simulations re¬ 
sults for the magnetization, sub-lattice order parameters, 
specihc heat and susceptibility are presented in Sec. IV. 
A summary, conclusions and directions for future work 
are given in Sec. V. 


II. THE MODEL 


The dimensionless dipole interaction Hamiltonian for 
Heisenberg model spins is given by: 


T~('dip — 
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where fij is the (dimensionless) vector connecting spins 
on lattice sites i and j. This Hamiltonian contains two 
terms that offer some insight into how the spins may tend 
to align at low temperatures. The second term in Eq. [T] 
favours spins that align along the vectors Therefore, 
in the ground state it is reasonable to expect spins to 
exhibit some preferred orientation related to the lattice 
vectors. This phenomenon is referred to shape anisotropy 
and is dependent on the geometry of the lattice. The 
first term is noticeably similar to the antiferromagnetic 
exchange interaction. 

In view of what is known about the ground state of the 
exchange-only kagoine lattice, the ground state of the 
dipolar kagome system can be expected to be a three- 
sublattice system with fixed spin orientations for each 
sublattice. We expect that in the ground state all spins 
on sites of the same sublattice will point in the same 
direction, and so the ground state can be represented by 
three angles, one for each sublattice. From the above 
arguments it can be seen that in the ground state, the 
spins will tend to lie in the plane of the lattice due to the 
shape anisotropy induced by the dipole interaction. In 
order to efficiently compute the long-range interaction, 
we use standard Ewald summation techniquesi^ 


III. THE GROUND STATE 

The EFM has been successful in the determina¬ 
tion of complicated spin structures associated with spin 
glassesi^ In our implementation of this method, the 
starting point is an ordered or random initial spin con¬ 
figuration, the spins of which are iteratively aligned with 
the local field to minimize the energy. This procedure is 
then repeated multiples times to ensure that the global 
minimum energy spin state is achieved. In addition, low 
temperature MC simulations presented in SectionlTVIcor- 
roborate the ground state results. 
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FIG. I: Example low-temperature spin structure iden¬ 
tified by MC and EFM simulations. The state contains 
states 1 & 5 as defined in Table ID Domain walls sepa¬ 
rate the two ground states. Spins in a domain wall have 
an orientation that is the average of the orientations of 
the spins belonging to the sublattice on which the do¬ 
main wall forms. The top two rows are state 1 while the 
bottom rows correspond to state 5. 
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(b) State 2 



(c) State 3 (d) State 4 
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(e) State 5 

(f) State 6 


FIG. 2: Six pure ground state domains of the 2D dipolar 
kagome spin lattice. 


This analysis reveals multi-domain spin configurations 
involving pairings of six types of magnetic states, sep¬ 
arated by domain walls as shown in Fig [T] These six 
ground states, defined in Table HI all exhibit a three- 
sublattice structure with fixed spin orientations for each 
sublattice. Each state exhibits shape anisotropy as pre¬ 
dicted, as one sublattice in each state is aligned along the 
direction oi 0 = ^ relative to the horizontal axis where 
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n is an integer, n G {0,1, 2, 3,4, 5}. The orientation of 
spins on the other two sublattices deviate from 9 by an 
angle (j) = 36.3887° such that the angles of the other two 
sublattices are given hy 9± = 9 zL (p, as seen in Table I. 
The energy of each of these six spin configurations is the 
same, giving rise to a six-fold degeneracy. The six ground 
state spin configurations are shown in Fig. [21 


rations of the states reported by Tomita and calculating 
the energy of the system, we find that the energy per 
spin, £i/spin Ri —1.6 to be greater than the energy per 
spin for our ground state, Fi/spin = —2.38895. We note 
that the transition temperature found by Tomita is es¬ 
sentially the same as ours (T^r ~ 0.43, as shown in the 
next section). 


Domain 

Oa 

Ob 

Oc 

1 

23.6113 

96.3887 

60 

2 

203.6114 

276.3887 

240 

3 

36.3887 

0 

-36.3887 

4 

216.3887 

180 

143.6114 

5 

120 

83.6113 

156.3887 

6 

300 

263.6113 

336.3887 


TABLE I: Ground states for classical dipoles on the 
kagome lattice. Note that the number of the state does 
not correspond to the value of n that defines the overall 
orientation of the magnetization of that state. 

To further illustrate how the ground states combine to 
form domain walls at low temperatures, Fig. [3] shows 
spin configurations 3 and 6 from Table U as pure ground 
states and then shows the mixture of these two states. 

The ground states have an overall orientation that can 
be defined by considering a “macrospin” formed by three 
spins located on a single triangle within the lattice. The 
overall orientation 9m of a domain is defined by the aver¬ 
age orientation of the three spins that form a macrospin, 
Gm = [9 + 9+ + 9-)/2> = 9. Therefore, each state can 
be defined by an orientation 9m = ^,n G 0,1, 2, 3,4, 5. 
Our shape anisotropy predictions made from the mag¬ 
netic dipole interaction term are confirmed further as the 
spins all lie in the plane of the lattice with deviations from 
the plane being minimal. 

At low temperatures, the system orders itself into mix¬ 
tures of these states separated by domain walls that have 
a very small effect on the systems total energy. There are 
six favourable pairs of states at low temperatures: (1, 5), 
(1,3), (2,6), (2,4), (3,6) and (4,5). Within these spe¬ 
cific state pairings, the two participating domains have a 
sublattice with similar orientation, and for each of these 
pairs, 9m in the two domains differs by 60°. For example, 
states 1 {9m = 60°) and 5 {9m = 120°) have similar val¬ 
ues of 9b, as seen in TablelH The triangular lattice, which 
orders ferromagnetically, has less geometrical frustration 
than the kagome lattice, which may be responsible for 
the difference in ground state spin structures. 

In a previous MC study on the dipolar kagome sys¬ 
tem, Tomita suggested that the spin structure of the sys¬ 
tem at low temperatures is composed of ferromagnetic 
chains with spins between the chains appearing to not 
be ordered^. In fact, snapshots of the spin-structure 
from our simulations at temperatures studied by Tomita, 
indicate distinctly different spin structures. As a result 
of this disagreement, we develop a simple computational 
model to verify our results. By generating many configu- 
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(b) State 6 
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FIG. 3: Diagram depicting domains 3 (black arrows) and 
6 (red arrows) and the domain wall (blue) formed by the 
mixture of these states. Sublattice sites are coloured blue 
(A), red {B) and green (C). The spins separating the two 
states form the domain wall with an orientation that is 
the average of the orientations of the spins on sublattice 
C (green sites) from states 3 and 6. 
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Because the magnetic dipole interaction is a long-range 
interaction, it is of interest to study how the ground state 
structure evolves as a function of the range of interaction. 
By considering a finite-size model of three spins belong¬ 
ing to a central cluster on a kagome lattice, as in Fig. |4l 
we can study the effect of the range of the interaction on 
the ground state of the system. For simplicity, we have 
assumed that the ground state is a three-sublattice sys¬ 
tem so that all spins belonging to a sublattice have the 
same orientation. Only interactions between the central 
spins and their neighbours are considered as shown in 
Fig. m That is, the spins outside of the central cluster 
do not interact with each other. The range of interaction 
is varied at each step and the energy is minimized nu¬ 
merically in three variables, 9a, Ob and 9c, the angles of 
spin orientation for sublattices A, B and C, respectively. 



FIG. 4: Visualization of the hnite-size model. Spins be¬ 
longing to a central cluster (denoted by red, blue and 
green lattice sites) interact only with their neighbours 
within the range of interaction. 



2 4 6 8 10 12 


FIG. 5: Sublattice angles as a function of the range of 
interaction. The angles slowly converge as the range of 
interaction is increased. The ground state is nearly real¬ 
ized after including up to only second nearest-neighbours. 


The plot of the sublattice angles as a function of the 
range of interaction shown in Fig. [S] indicates that the 
underlying physics of the magnetic dipole interaction is 


captured after the range includes the second nearest- 
neighbours of the central spins. Interestingly, when the 
range only includes the nearest-neighbours we retrieve 
a result that is similar to the q = 0 ground state of 
the exchange-only kagome latticeJ^ That is, the nearest- 
neighbour dipolar interaction gives a 120° spin structure 
with the total magnetization on a triangle equal to zero 
coupled with the shape anisotropy as all the angles point 
along lattice vectors. Note that in this model, the spins 
in the central cluster interact with each other twice. This 
double counting reflects the Ewald summations used in 
our simulations to calculate the correct energy per par¬ 
ticle of the infinite periodic system. 


IV. SIMULATION RESULTS 


Results from using Metropolis MG simulations at fi¬ 
nite temperatures are described here. All quantities cal¬ 
culated are averaged over 10® MG steps (with the initial 
10% discarded for thermalization) unless otherwise spec¬ 
ified and are calculated for different lattice sizes L x L 
with N = 3L^/4 spins. Here, L represents the number 
of points along an edge of the underlying triangular lat¬ 
tice, with one-quarter of the sites removed. A single MG 
step refers to N individual Metropolis spin update at¬ 
tempts. Simulations versus temperature are done three 
different ways>^ as cooling runs with an initial random 
spin configuration and then using the final configuration 
of the previous temperature as the initial configuration 
of the next lower value of T ; heating runs with the ini¬ 
tial configuration being one of the six ground states; and 
finally independent temperature runs where a random 
initial configuration is used at each value of T. Results 
for the energy calculated from cooling simulations with 
AT « 0.01 are shown in Fig. [Slfor L = 18 and 24 showing 
little impact of system size. There is a small inflection 
in the energy at T Ri 0.4, indicating the possibility of a 
phase transition. 



FIG. 6: Energy as a function of temperature for different 
lattice sizes L from MG cooling simulations. 
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To determine if the inflection in the energy corresponds 
to a phase transition, the specific heat (per spin) as a 
function of the temperature is calculated from cooling 
runs. Fig. 0 shows a peak at T fv 0.43, consistent with 
Tomitai^^, with a relatively weak dependence on system 
size that may be indicative of a continuous phase transi¬ 
tion. 
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FIG. 7: Specific heat versus temperature for different 
lattice sizes L from MC cooling simulations. 


In classical spin systems, the total ferromagnetic mag¬ 
netization, 


Mf = 



( 2 ) 


is a useful quantity to consider as it can offer insight 
into the order of the system. It Mf =0, the system is 
completely disordered with no net magnetization, while 
if = 1 the system is in a ferromagnetic phase with 
all the spins aligned along the same orientation. Fig. [5] 
shows the total ferromagnetic magnetization from cool¬ 
ing runs and indicates that the system is in a disordered 
phase at temperatures above Tpf. At low temperatures, 
the results suggest that the system is in a state that ex¬ 
hibits order with some deviations from fully ferromag¬ 
netic order. This is expected from the analysis presented 
in Section uni as low-temperature states are composed of 
a mixture of ground states which are not totally ferro¬ 
magnetic, that is Mf ^ I for a pure ground state. In 
fact, for a pure ground state, Mf = 0.8700, which is 
greater than the value found for an equal mixture of two 
domains, for which Mf = 0.8700 cos 30° = 0.7534 (since 
the macrospin angles in the two domains are separated 
by 60°, as discussed in Section ED). 
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FIG. 8 : Total ferromagnetic magnetization M/ as a func¬ 
tion of temperature for various lattice sizes from cooling 
runs. 



Based on the results in Section lllll we expect that each 
sublattice is fully ordered in the ground state. Therefore, 
we calculate the sublattice magnetization order parame¬ 
ter. 


M^ = 



( 3 ) 


where 7 represents sublattice A, B or C, and k runs 
over all 7V/3 spins of the sublattice; as well as the total 
sublattice magnetization order parameter, 



( 4 ) 


When the value of this total sublattice magnetization is 
one, the system has three ferromagnetically ordered sub¬ 
lattices as seen in the ground states defined in Section ITTll 
However, at low temperatures the presence of mixtures 
of ground states separated by domain walls will reduce 
the value of Mt from unity to a value that depends on 
the number of spins that belong to each domain present. 
Fig. ini confirms this prediction as the total sublattice 
magnetization is not equal to unity at very low tem¬ 
peratures. For a state consisting of an equal number of 
spins belonging to each domain, Mt = 0.7757. In Fig. [HI 
Mt tends to approximately 0.78, which is slightly higher 
than if there were and equal number of spins belonging to 
each domain. This is expected as the domain with more 
spins will contribute more to the magnetization, bringing 
the value closer to the pure ground state value of unity. 
Owing to the small energetic penalty of forming domain 
walls, the ground state is not realized as the appearance 
of a second domain will occur even at the lowest of T. 
The energy difference between a pure ground state of size 
L = 12 and one with a single row of a complementary do¬ 
main is AE = 0.38 (AA/spin r; 4 x 10“^ and AE scales 
linearly with L). This excess energy is so small that mix¬ 
tures of states are inevitable in an ergodic system. By 
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contrast, a single row of a non-compleinentary domain 
(e.g. state 2 with a row of state 1) yields an energetic 
penalty two orders of magnitude higher. 


0.7 
0.6 
^^0.5 
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T 

FIG. 9: Total sublattice magnetization as a function from 
cooling runs. Small fluctuations at T ~ 0.3 correspond 
to the system “deciding” which combination of domains 
will be present at low T. 



X / \ \ / X \ » X \ \ / \ X X 

X / t i \ — i 

X / X i X XX \ x\>.XXX \ \ X — 

/ X / \ / X j „ 

X X \ X X X X — — ( X X X X X X 

— / - J \ i i ^ 

X / \ X X X X X X X \ X XX \ r 

X / X \ X / .. 

\ / X X X X \ X \ ! X . X X \ • X / 

X / X \ \ \ , t - 

X / \ X X X \ X X \ X X X X \ X X 

X t i X \ -X / — 

xx-'-^x XXX, Xxxxx-xxX^j XX 

X X i/x X X Fx X \ \ \ X \ \ \ X 
t -- \ / i 'X 

--\ X / X \ X ^ \ X xx \ X / 

- / ■- i \ ( - \ - 




'^Y 

.y 


(a) T = 0.3176 


In both Fig. [5] and Fig. O there exists a common fea¬ 
ture at r ft! 0.27 where the magnetization parameters 
appear to fluctuate. Because these fluctuations are only 
present for a narrow region of temperatures, we examine 
the possibility that they are associated with an onset of 
a reduction in thermal-fluctation induced switching be¬ 
tween degenerate spin states. We note that there is no 
discernable feature at this temperature region in the spe¬ 
cific heat. 

Fig. [TUI shows that the sublattice ferromagnetic mag¬ 
netizations for each sublattice of the system have large 
fluctuations in the region where the features are present 
in Figs. |H]and[Ul 



FIG. 10; Total ferromagnetic magnetization for indi¬ 
vidual sublattices A, B, and G from cooling runs with 
L = 18. Large fluctuations correspond to the sys¬ 
tem changing state mixtures until the system becomes 
“frozen” into a combination of states. 


X \ X \ X \ \ \ \ X \ X \ X \ X \ 

\ 1 / \ f I ^ X * 

/ / X t X / X X / X / X X X / .. 


I 


\ ! f / I 


-Axxxy x\ \ \'xxx\xxx \ 
t f 1 ’ f ^ 1 / t 

~- \ X — \ X \ \ \ x X X \ ~- \ X \ 

/ t / t f ‘ \ / 

X - \ / fx/-xx/x.--.tx/.-- 

f f‘ f \ \ \ ’ ’ 

X \ X \ \ \ t XXX-^ \xx \x \ 

\ \ ^ f \ t ^ * t 

/X • X x X / X / -— / X / X / — / X 
^ T \ 1 ^ r / \ ^ 

! X/X/ xX—/ —XX/X—X/X 

t f / t * / / f 

X\— \ XXx\xXx\x > x\— ^ 

i f \ i f / / i / 

YYYYYYYYY 
YYYYYYYYY 
YYYYYYYYY 
YYYYYYYYY 
YYYYYYYYY 
YYYY YY'^YY 
YYYYYYYYY 

Y-y YYYYYYy 

(b) T = 0.2067 


FIG. 11: Snapshots of the spin structure in the region 
where fluctuations in the magnetization order parameters 
are present (See Fig. [TOl) . Below the snapshots, we have 
included diagrams with coloured triangles used to easily 
identify which states are present. Each state has an asso¬ 
ciated colour: State 1 (Magenta), State 2 (Yellow), State 
3 (Red), State 4 (Green), State 5 (Blue), State 6 (Or¬ 
ange). The colour of each triangle is calculated using the 
orientation of a macrospin Bm as a variable in a continu¬ 
ous “colour function” that associates a Red/Green/Blue 
value with Bm- a) States 3 and 6 are present, b) States 
1 and 5 are present. 
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By examining the spin structure in this temperature 
region (Fig. [TT|) the large fluctuations can be explained 
by the system changing which state pairs are present. In 
Fig. [TTh . a snapshot of the MC results at a higher tem¬ 
perature, T = 0.3176, shows states 3 and 6 are present 
while at T = 0.2067, shown in Fig fTTb . states 1 and 5 
are present. For T ^ 0.3, the energy required to change 
states is low enough to allow for large fluctuations be¬ 
tween spin configurations. Once the system is cooled 
to a sufficiently low temperature, the energy required to 
change from one pair of states to a different pair becomes 
large enough that the probability of the system changing 
states becomes very small. At this temperature T ^ 0.2, 
the system becomes “frozen” into a configuration that 
consists of a particular pair of states. 


That two sublattice magnetizations are significantly 
less than unity and one sublattice magnetization is nearly 
at unity (as in Fig. uni) can be understood by considering 
Fig. [3l where spins on sublattice C all nearly point in 
the same direction, while for sublattice A (and B) the 
spins are offset by 96.3887° in the two domains, giving 
Ma = Mb = cos (96.3887/2) = 0.6666. 


In Fig. [T21 results are presented for cooling, heating 
and individual temperature simulations showing the total 
ferromagnetic magnetization as a function of tempera¬ 
ture. Because heating simulations start from a ground 
state, the ferromagnetic magnetization is higher than 
that of a mixture of domains which have spins with op¬ 
posing magnetizations. As the temperature is increased, 
there are fluctuations introduced to the ground state spin 
orientations that lower the magnetization. At T Ri 0.2, 
the magnetization sharply drops to values consistent with 
the cooling simulation results. Similar to the system 
becoming “frozen” into a mixture of domains when it 
is cooled to T ^ 0.2, the system “melts” into a state 
composed of a mixture of domains when it is heated to 
Ta; 0.2. 


Also shown in Fig. [T^are results from individual tem¬ 
perature simulations which have been averaged over 10^ 
MC Steps. The magnetization is in agreement with both 
cooling and heating simulations for temperatures T > 0.2 
but has a large amount of noise for temperatures below 
T ft! 0.2. This noise is resultant from the fact that at 
each temperature studied the system has to be equili¬ 
brated from a randomly generated initial state. Since 
the acceptance probability is proportional to the Boltz¬ 
mann factor, as the temperature is decreased the accep¬ 
tance probability will decrease exponentially. Therefore, 
at low temperatures the number of Monte Carlo Steps re¬ 
quired to equilibrate the system increases greatly which 
increases the computational time required to perform a 
simulation that yields good statistics. 
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FIG. 12: Total ferromagnetic magnetization from cool¬ 
ing, heating and individual temperature simulations on 
lattice of size L = 18. 

Finally, Fig. IT^ shows the magnetic susceptibilityi^ as 
a function of temperature for cooling and heating sim¬ 
ulations. Both simulations yield a peak at T ft 0.43 as 
expected from previous results that identified a phase 
transition at this temperature. The results also identify 
features corresponding to how the system changes as it is 
heated or cooled. From the heating results, we see a small 
peak in the susceptibility at T ft 0.2 where the system 
begins to form domain walls. Cooling simulation results 
show that the susceptibility experiences fluctuations at a 
temperature range close to T ft 0.3 which results from 
the system reordering into mixtures of domains during 
this range. 
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FIG. 13: Magnetic Susceptibility as a function of tem¬ 
perature for cooling and heating simulations on a lattice 
of size L = 18. 
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V. SUMMARY AND CONCLUSIONS 

As a prelude to studies of magnetic thin-film ge¬ 
ometries which mimic spin-valve exchange bias struc¬ 
tures, and as a compliment to studies of magnetic nano- 






structures, classical Heisenberg spins with only dipole- 
dipole interactions on the 2D kagome lattice have been 
studied to identify the ground state spin structures, de¬ 
generacies and thermal behaviour of this geometrically 
frustrated lattice. In contrast with the simple triangu¬ 
lar lattice with edge-sharing triangles which exhibits a 
ferromagnetic ground state, our EFM and low-T MC 
simulations reveal that the corner-sharing kagome struc¬ 
ture leads to a more complicated six-fold degenerate spin 
states. These have been characterized (see Table I) in 
terms of three ferromagnetically ordered sublattices and 
multi-state domain structures that result from the simu¬ 
lations have been identified. Similarities and differences 
from the near-neighbor antiferromagnetic exchange only 
120 ° spin structures of the triangular and kagome lat¬ 
tices are described and result from beyond-third-neighbor 
dipolar terms where shape anisotropy effects become 
more important. 

The thermal behaviour of the system was studied 
through cooling, heating and individual temperature MC 
simulations. Results suggest that the system undergoes 
a phase transition at T Ri 0.43 in agreement with pre¬ 
vious MC simulations^^ but the nature of the ordered 
state differs. Ferromagnetic and sublattice magnetiza¬ 
tion order parameters calculated during cooling simula¬ 
tions suggest that below the critical temperature, the sys¬ 
tem changes between mixtures of ground state domains 
with fluctuations present until a threshold temperature 


T Ri 0.2 is reached where the system becomes “frozen” 
into a specific mixture of ground state domains. Heating 
simulations yield corroborating results in that the system 
has only one domain present at low temperatures until it 
reaches T r; 0.2 where the system readily changes into a 
mixture of domains as reflected in Fig. [T^ Our results 
suggest that the phase trasnition at Tjv is continuous and 
its criticality is likely impacted by domain fluctuations 
and the long-range nature of the dipole interaction^ 

Also of interest is to extend this analysis to study 
dipole interactions on the 3D fee ABC-stacked kagome 
systemsJ^ Preliminary results^E indicate a phase transi¬ 
tion to long range order at ~ 0.4 to a state with a 
different spin structure from the 2D kagome system and 
also different from the ferromagnetic order in the stan¬ 
dard fee lattice 

Finally we note that during the preparation of this 
manuscript we became aware of work that shows similar 
six-fold degenerate ground states as described here^S 
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